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Abstract 

We show theoretically that the emission spectrum of a single large quantum dot strongly coupled 
to a single photon mode in a microcavity can be qualitatively different from the spectrum obtained 
with an atom in a cavity. Instead of the well-known Mollow triplet we predict appearance of 
multiplets with the number of peaks a function of the quantum dot size and pumping intensity. 
The mutiplets can appear if the quantum dot is larger than the exciton Bohr radius, so that 
excitons are confined as whole particles in the dot. In this case the Pauli principle is relaxed and 
one can accommodate more than one exciton (but still a finite number) in a given quantum state. 

PACS numbers: 42.50.Ct, 78.67.Hc, 42.50.Pq 
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Introduction — Coupling of a zero- dimensional photonic mode to a zero-dimensional ex- 
citon (electron-hole]state has been realized in pillar and photonic microcavities containing 
quantum dots DDDD . Very recently, two groups have obtained independently the experi- 
mental evidence for vacuum-field Rabi splitting in such structures |3j, \$ . These observations 
confirmed theoretical predictions of strong exciton-photon coupling in zero-dimensional sys- 
tems BQ. 

Nonlinear optical effects in strongly coupled zero-dimensional systems are expected to 
be extremely rich. For a while, only a manifestation of the Purcell effect hjl had been 
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experimentally reported |2|, while theoretical predictions went much beyond, evoking 
appearance of dressed exciton and bi-excitons states sj]. The crucial question for description 
of emission from quantum dot cavities (QDC) is whether crystal excitations coupled to 
light behave like fermions or like bosons in these systems The fermionic regime has 
been a subject of cavity quantum electrodynamics (QED) for decades. It is familiar in 
atomic cavities, and can be efficiently described by Dicke formalism Q]. Its signature is the 
appearance of the Mollow triplet [jjj in emission spectra of the cavity. 

Clearly, the fermionic regime is realized if carriers are strongly confined in small-size 
quantum dots (case of for example). On the other hand, it does not hold if the dot size is 
much larger than the exciton Bohr radius (further referred to as the "large dot" case, realized 



in particular in the dots originated from islands of quantum well width fluctuation 0, 
In large dots, excitons are quantized as whole particles and exhibit bosonic properties at 
least to some extent. The limiting case of a large dot is a quantum well where quasi-Bose 
condensation of excitons is indeed possible 1131 and Rabi-doublets are observed in the spectra 

n 

of microcavities 114. Dicke formalism does not hold in this it does not take into 

account the fundamental difference between an atomic and an excitonic system: the number 
of atoms is fixed while the number of excitons is not. From a formal point of view in the 
purely fermionic case the light mode is coupled to a two-level system, while in the purely 
bosonic case it is coupled to a harmonic oscillator having an infinite number of equidistant 
energy levels. Logically, the straightforward intermediate case would correspond to the 
optical mode coupled to a truncated harmonic oscillator, having a finite number, say N, of 
energy levels. As such, a large QD offers a realization of parastatistics which interpolate 
between Bose-Eistein and Fermi-Dirac statistics by allowing up to N > 1 particles in a given 
quantum state jlfl llfil]. 
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In this work we study theoretically the emission of light by a QDC using the model of a 
truncated harmonic oscillator. We assume that a finite number of excitons (N > 1) can be 
generated in the quantum state of interest. Here N is taken as a free parameter, expected to 
scale like S^/ita^ where S^ ot is the cross-section surface of the dot and a-Q the exciton Bohr- 
radius. Physically, it corresponds to the following image: the dot is much larger than the 
exciton Bohr radius, so that excitons are confined as whole particles inside it. N = 1 recovers 
the fermionic case. We also assume that the energy separation between quantum confined 
excitonic states exceeds the vacuum Rabi splitting, so that one can neglect all excitonic 
states except one, strongly coupled to a single light mode. This latter condition, though very 
artificial for quantum well microcavities, is easily met in QDC, where typical values of the 

, we neglect exciton-exciton interactions, 
which allows us to assume the energy of the excitonic transition to be independent on the 
number m of excitons already created in a given quantum state if m < N. If m = N 
no more excitons can be created in the quantum state we consider. Finally, we assume 
the exciton-photon system to be fully spin-polarized, so that polarization effects can be 
neglected and bi-excitons bound state impossible. Experimentally this corresponds to the 
resonant circularly polarised pumping and no efficient spin-relaxation time in the system. 
In general to deal with QDC one should solve the complicated many-particles problem 
accounting for exciton-exciton, exciton-carrier and exciton-photon interactions. Here we 
deal with the simplified model described above as an ideal limit of more refined microscopic 
calculations which display quantitative departures but an overall qualitative agreement. 
Our present model retains the essential physical difference of a large quantum dot from a 
small one and a quantum well and allows to obtain analytically the emission spectrum of 
the structure. We show that this spectrum is remarkably different from well-known Mollow- 
triplet and Rabi-doublet. It exhibits multiplets whose number is sensitive to the dot size 
and pumping intensity. 

Formalism — The electromagnetic field is modelled in the approximation of the quasi- 
mode coupling by single mode Bose field annihilation operator a. We introduce for a large 
quantum dot, able to host N excitons, the operator o~n which behaves like a Bose field 
annihilation operator when the number of excitations is less than or equal to N, i.e., 

(Jtv \m) = y/m \m — 1) , crjy \m) = \/m + 1 \m + 1) m ,jv 
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where 9 m ^ is the modified Heaviside step function which is zero if m > N and one if m < N. 
Its matrix representation in the basis of Fock states is the N + 1 squared matrix which is a 
truncation of the annihilation matrix of a Bose field. For instance, for N = 2, 

^0 1 ^ 

o 2 = V2 ■ (1) 
^0 j 

The algebra for this operator follows straightforwardly: 

[o~n, cr N ] = ljv+i — Atv+i , (2) 

where Ijv+i is the N + 1 squared identity matrix and (Ajv+i)^- — (N + 1)6^6^ (with < 
iJ<N). 

The hamiltonian for our system is similar to Dicke [l0|] or Jaynes-Cummings [l7| hamil- 
tonian but instead of coupling the field a with a spinor J or a two-level system (fermion) a\, 
one couples a with this operator a N , also in the rotating wave approximation: 

H = hjja)a + hjja^ N (jN (3a) 
+ hg(a)<y N + aa^ N ) . (3b) 

We considered zero detuning for simplicity The general case with a small detuning brings 
fine quantitative departures which are the topic for a separate communication. Here hp is the 
coupling constant corresponding to (half) the Rabi splitting in the linear regime ^, |j|. First 
line ()3a)l describes bare states of photons and excitons, the later being limited to a maximum 
N. We note eigenstates of this first line alone \n, m) with n the number of photons, any 
integer number, and m the number of excitons, any integer less than or equal to N. Second 
line (|3b|) describes the coupling between the two fields which in absence of dissipation will 
always be "strong" in the quantum mechanical sense. The dissipation is caused by two 
channels of energy escape included as non-hermitian contributions of the Lindblad type to 
Von Neumann equation, 

ihp = [H, p] + L a p + C aN p , (4) 

where 

L A = - y(AUp - 2ApA* + pA^A) (5) 
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for A = a (resp. <Jn) describing the photon leakage through the cavity mirror by tunnel 
effect (resp. the emission from the active state within the cavity) with associated decay 
constant 7 a (resp. 7 CTJV ). The procedure is straightforward in principle: since the number of 
excitations is conserved by H, one diagonalises it in the subspace of n excitations and thus 
obtains new eigenstates of the system (so-called dressed states in cavity QED terminology). 
Thus, in the basis 7i n = {|0, n) , |1, n — 1) , . . . , \N, n — N)} (for n > N), the hamiltonian 
reduces to: 



+ y/(i + l)(n- 



(6) 



J 0<i,j<N 

where the term between square brackets is the generic expression for the interaction matrix 
(which is zero everywhere but below and above its diagonal). Case N = 1 gives the Mollow 
triplet. 

Along with general results, we will deal explicitely with the all-important case N = 2. 
It is expected the effect we predict will first be observed for this very value N = 2. In this 
case, the matrix representation (JHJ) in the basis H.2 reads for n excitations (with n > 2) 
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H = ftwnlz + hg 
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(7) 



yO y^n-V) U j 
Cases where n is lower than or equal to N fall into linear regime displaying two peaks, the 
exciton field never reaching its saturation density. If n > N, the new N + 1 eigenstates 
(dressed states) will be denoted, for the manifold n, by \u) with v = 1, . . . , N + 1 indexing 
the state. For N = 2, 
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\» = 3> B = 

in the bare states basis 7^2- The associated eigenvalues are: 

El = hum + (u- 2)hgV3n - 



(8a) 
(8b) 
(8c) 

(9) 



The spectrum of emission comes from transitions between manifolds with n and n — 1 
excitations. One can obtain the exact expression for the luminescence spectrum from (J3J) and 
quantum regression theorem to evaluate the time correlation function needed for Wiener- 
Khintchin theorem. We content here with Fermi's golden rule formula for perturbations 
responsible for (jSJ). These perturbations arise from the weak coupling of cavity photons a 
(resp. excitons <tjv) with the continuum of external modes of the electromagnetic field in 
vacuum state. Therefore, calling c"' n the projection of state \ct) n on state \i,n — i) one can 
obtain the coefficients in (jHJ) as: 

N 
i=0 

The corresponding transition probabilities between states \u, n) and |i/,n— 1) are obtained 
as: 

N 2 

E, u'n—l\* v,n r 
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Mn-l^|f) n | 
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In cavity QED terminology, 7 en d and Ji at correspond to end-emission and lateral-emission 
photo-detection respectively, while in luminescence of a microcavity, one observes the sum 
of the two contributions 7i at + 7 e nd simultaneously. To detect I\ at independently one should 
use the scattering geometry. 

Results — Our main result is the appearance of multiplets in the configuration n > N. 
This is in stark contrast with cavity QED where nonlinearity results in apparition of the 
Mollow triplet [ll]]). Three questions need be solved in the general case: first the number of 
peaks allowed by the transitions which is of a combinatorial nature, providing the number of 
transitions between two multiplets taking into account their degeneracy. Second the energy 
splitting between the various peaks thus obtained which is the eigenvalues problem. Third 
the strength of transitions which is the eigenvectors problem. If the radiative broadening is 
so small that all peaks can be discriminated and if no selection rules forbid some transitions, 
the multiplets structure is that of 2N + 1 groups of more or less closely packed peaks for 
a total number of (N + l) 2 . We label fi = the group at the centre of the spectrum, 
and /i = ±1, ±2, . . . groups located symmetrically about the central one. The number of 
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peaks is N + 1 — \fi\ in the /i th group. Each peak is identified unambiguously by a couple 
of integers (k, I) with < k,l < N. Peak (k, I) belongs to group \i — k — I and is the k th 
(resp. Z th ) if located at the left (resp. right) of the spectrum, counting as first in each group 
that of lower energy. The energy of this peak E(k, I) is given by the photon energy released 
in the transition between two manifolds and is therefore given by the difference of initial 
and final eigenvalues E% in the dressed states basis: 

E(k,l) = El - E?' 1 (12) 

In the case N = 2, from Q and (|12|) we get in good approximation (the better the higher n) 
the mean energy of peak (k, I) as 

E{k, I) « Huj + {k - l)hgV3^ - (2k -51 + 6)-^L (13) 

2v3n 



This gives the structure of the multiplet: there is a splitting of magnitude \u\hg\/3n between 
the central (fi = 0) and the i4h group of peaks. The central group has a fine structure 
splitting (between (1,1) and (3,3)) of \pihgj ' y/n. At the onset of the nonlinear effect (n = 3) 



the splitting is close to the Rabi splitting 2yJ(hg) 2 — — r ) cr2 ) 2 The fine structure 
splitting of the groups \x = ±2 is half the central one: V / 3%/(2 v / n). As n increases, each 
group of peaks gets farther from the others while inside this group, the fine structure splitting 
is reduced. This behaviour is seen in Fig. (0). When the photon field is so intense that for 
the given broadening 7 CT2 the fine splitting is unobservable, the spectra resemble a Mollow 
triplet. It can be distinguished from it by the splitting which is v^3/2 times smaller and the 
presence of additional (small) peaks with a splitting \/3 larger. Also the ratio between the 
central peak and its closest satellites is 9/4 rather than 1/2 in case of the Mollow triplet 
(see Fig. flUo)). 

This structure applies both for 7 en d and Ji at spectra. However, they display two drastically 
different behaviours. First J cn d is maximum at the center of the spectrum with peak (2,2) 
the highest one. The situation is opposite with Ji at for which, whatever n, peak (2,2) is an 
exactly forbidden transition. 

The same procedure can be applied to larger values of N, e.g., for N = 3 (cf. |HJ): 



= hum + sign(z/ - -)hg\J (-1)" 'v^V '3n 2 - 9n + 8 + 3n + 4 (14) 
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FIG. 1: Multiplets observed in the spectrum 7i at of lateral emission of the cavity for (a) N = 2 
and n = 3 and (b) n = 5 (solid) and n = 10 (dashed) (also for N = 2). The coupling constant hg = 
50/xeV. Broadening is obtained by convolution with Lorentz functions of width j a = 100/xeV 
and 7 CT2 = 5/ieV |3j. The inset in Fig. (a) showns J e nd f° r n = 3, which is a Lorentz function; the 
inset of Fig. (b) shows ii a t in the limit of huge n (n = 100). 

with 1 < v < 4, and for iV = 4: 

.E™ = hnu + sign(z/ — 3)hgx 

(15) 

sign(3 - v)(-iyV2Vbn 2 - 25n + 38 + 5n - 10, 

this time with 1 < v < 5, etc... (where sign(x) = if x = and is x/\x\ other- 
wise.) Along the same lines as before, one can derive the structure and splittings of the 
spectrum of emission in these cases. For instance, for N = 3, the (l,l)-(3,3) splitting 
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FIG. 2: Multiplets for (a) N = 3, n = 5 (inset N = 4, n = 5) and (b) iV = 20, n = 25 with same 
parameters as in Fig. 

is (v3 — V6)hg j 1 \/n rs 0.74hg/ y/n. Note that the difference in energy enters as square 
roots inside a square root and thus becomes smaller and smaller with increasing N. Fig- 
ure (J2t>) shows the case of an extremely large dot where satellite peaks become suppressed 
and the Rabi-like doublet dominates the spectrum. For N — > oo the familiar Rabi-doublet 
seen in the spectra of quantum well cavities is recovered by this model. 

Conclusions - - We predicted appearance of multiplets in emission spectra of zero- 
dimensional photonic cavities containing single large quantum dots. The effect comes from 
relaxation of the Pauli principle for excitons confined as whole particles in large dots. The 
effect is strictly nonlinear as it requires the number of photons injected in the cavity to 
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exceed the capacity of the dot. This capacity governs the number of peaks in emission. In 
the small dot limit, it reduces to three (Mollow triplets). 
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